Vortex lattice in presence of a tunable periodic pinning potential 
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The vortex patterns stabilized by the square array of artifi- 
cial pinning sites with a tunable pinning strength are studied 
by using a phenomenological approach in the London limit. 
The transitions between pinned and deformed triangular lat- 
tices are analyzed as a function of the amplitude of the vortex- 
pinning site interaction and the characteristic length-scale of 
this interaction. The critical current and different phases of 
vortex lattice are studied in presence of external transport 
current. 

PACS numbers: 74.20.De, 74.60.Ec, 74.25.Ha 



I. INTRODUCTION 



over the pinning forces and the triangular lattice can be 
recovered. 

In the present paper we have focused on pinning phe- 
nomena in presence of a square array with a tunable pin- 
ning strength. We reveal how stable vortex configura- 
tions are evolving as a function of the pinning strength. 
We have found the transition lines between pinned and 
deformed triangular vortex lattices. At certain matching 
fields the vortex lattice phase diagram is more compli- 
cated due to the existence of specific periodic phase. We 
analyzed the phase diagrams taking into account differ- 
ent types of pinning potential profile. The critical cur- 
rent in this phase is calculated and the phase diagrams 
are obtained. 



II. PHASE DIAGRAM 
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Nanoengineered periodic pinning arrays (PPA) have 
attracted recently a lot of attention [0-|9j. By using 
modern e-beam lithography, very well defined periodic 
arrays of submicron antidotes |j],|2|,||-^| or magnetic dots 
can be fabricated in a superconducting film with 
a low intrinsic pinning. The latter is needed to have in 
the film only artificial periodic pinning centers, which 
are well controlled. These laterally nanostructured films 
demonstrate sharp matching peaks in magnetization vs. 
field curves, corresponding to some specific vortex pat- 
terns stabilized by nanoengineered PPA. At tempera- 
tures quite close to the superconducting critical temper- 
ature, T c , the artificial pinning sites attract the vortices 
very strongly and, as a result, the underlying symme- 
try of the pinning array dictates the symmetry of sta- 
ble vortex patterns. For example, at the matching fields 
H = Hi and H — H\/2 (where Hi is the first matching 
field, which generates exactly one flux quantum per pin- 
ning site) , square pinning array stabilizes a square vortex 
lattice, while in a homogeneous non-patterned reference 
film a conventional triangular vortex lattice minimizes 
the vortex-vortex repulsive interaction. Stable vortex 
patterns in presence of the square pinning arrays were 
directly visualized with a help of the Lorentz microscopy 

If now the vortex pinning strength in a periodic square 
array is smoothly reduced, we can expect that sooner or 
later the vortex-vortex repulsion will start to dominate 



Let us consider the two-dimensional sample with a reg- 
ular square lattice of pinning sites with a period a. We 
assume here that the external magnetic field does not ex- 
ceed the first matching field Hi and is equal to one of the 
sub-matching fields. It is well-known that the regular tri- 
angular vortex lattice has the lowest energy in absence of 
any pinning ]Io| . Square lattice of pinning sites can im- 
pose its own symmetry on the vortex lattice. Therefore, 
vortex lattice can exist in the form of at least two sta- 
ble phases - pinned regular phase and triangular phase in 
two limiting cases (strong and weak pinning). The sym- 
metry of the pinned lattice is not necessarily the square 
one. Vortex patterns corresponding to different matching 
fields were listed in Ref. || for the limit of strong pinning. 
Some of these patterns we also present here on Fig. 1. 
When the pinning strength is small it is energetically fa- 
vorable to form a deformed triangular lattice. This lattice 
is nonregular, since triangular vortex lattice and square 
pinning site arrays are incommensurate. Transitions be- 
tween these phases were first studied in Ref. [|ll| in the 
limit of infinitely small length scale of pinning potential 
well. Here we will analyze general case of an arbitrary 
potential well size. 

The pinning energy density can be written phenomeno- 
logically as a Fourier expansion: 

... . \ -> 2-kvtlx 2imy , . 

U(x,y)=y A mn cos cos . (I) 
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In the pinned regular phase all vortices are located on 
the pinning sites and the energy of the system E\ (per 
one vortex) reduces to 



E 1 = E lat + U(0,0)=Ei at 



A 

r, 



(2) 



where U (0, 0) is the pinning potential in the center of the 
pinning site and Ei at is the energy of vortex interaction, 
which depends on the symmetry of the flux-line lattice 
(Fig. 1). 

The energy of the deformed triangular lattice E2 can 
be found in the framework of the elasticity theory, in 
which a discrete vortex lattice is replaced by an elastic 
medium. The accuracy of such an approximation will be 
discussed below. The energy E2 (per one vortex) is given 
by 



E2 = E tr 



E 



pin 



E, 



elast 1 



(3) 



where E tr is the energy of an ideal triangular flux-line 
lattice, Epi n is the pinning energy, and E e i ast is the en- 
ergy related to the lattice distortion. To find E pin and 
Eeiast we should calculate first the deviations of vortices 
from their positions in triangular lattice. We introduce 
two-dimensional vector of deviation of the elastic medium 
u(x,y). For E pm and E ela st we have 



Eeiast — 2g ' dS 



(Cll - C 66 ) (Vu) 2 + C 66 (VaU/3) 2 



Epi n = ^ J dS [U(x, y) - u(x, y)f (x, y)] 



(4) 



(5) 



where the integration is performed over the area S of the 
sample, Sq is the area of the unit cell of vortex lattice, 
and f = ~dll/dr is the pinning force acting on a vortex. 
In Eq. (5) we took into account only the first term in the 
energy expansion with respect to u(x,y). 

The deformation u(x, y) can be found by the minimiza- 
tion of the total energy (3). Using the Fourier represen- 
tation it is easy to obtain: 



Uka 



(kf k )k a k- 2 (kf fc )fc a - f k 

cnSok 4 



cnS'ok 2 



(6) 



where 



u k = J u(r)e lkr dS, f k = J f (r)e ikr dS*. 
Taking into account Eq. (1) one can find the following 
expressions for u(x,y): 



2ttchS{ 
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E^mn ~ 
1 



2nmx 2imy 
■ sin cos , (7) 



2irmx 2imy 
■ cos sin . (8) 



These expressions give the deviations of the vortices from 
their equilibrium positions in an ideal lattice. It is should 
be noted that the term proportional to 1/C66 in Eq. (6) 
disappears in Eqs. (7) and (8). This term corresponds 
to twisting of the vortex lattice. Absence of this con- 
tribution in the resulting expressions (7) and (8) implies 
that regular pinning potential does not twist the vortex 
lattice. 

Note that within the elasticity theory deviations must 
be much smaller than the intervortex distance d, which 
leads to the following constraint: 



< d. 



(9) 



Now we can find E e i as t and E p i n from Eqs. (4) and 



(5): 



1 



2cnS 



(< U 2 > - < U > 2 ) , (10) 



Eeiast =<U> (< U 2 > - < U > 2 ) . (11) 



Here we took into account the Parseval's identity 

Y / A 2 nn + A 2 = 2<U 2 > . (12) 



The elasticity energy defined by Eq. (10) is proportional 
to the dispersion of the pinning potential, which is a mea- 
sure of deviation of function from its average value and 
describes its sharpness. The first term in the right-hand 
side of Eq. (11) corresponds to the pinning energy in 
the zero approximation with respect to u(r). In this ap- 
proximation the vortex positions are independent of the 
pinning potential and the elastic energy equals zero. The 
last term describes the energy of the correlation between 
the vortex positions and the pinning sites. This term is 
equal exactly to —E e i a st and it is cancelled in the result- 
ing expression for the total energy: 



E2 — Et 



<U> 



(13) 



Note that Eqs. (11)-(13) are valid for any profile of the 
pinning potential. For illustration of our method we show 
the structure of a deformed triangular lattice calculated 
from Eqs. (7) and (8) in Fig. 2 at H u c n S /U o = 0.2, 
a = 0.1a. Open circles denote pinning sites, black dots - 
vortices. 

Comparing the energy of the pinned lattice E\ with 
the energy of the deformed triangular lattice E2 we can 
obtain the criteria for the transition between them. We 
choose the Gaussian profile for the vortex-pinning site 
interaction: 



U{x,y) = -C/ cxp -(r/a) 



(14) 



where Uo is amplitude of interaction, r is distance be- 
tween vortex and the pinning site, a and is the potential 



2 



length scale. The total pinning potential in the sample is 
equal to the sum of contributions from all pinning sites: 



U(x,y) = 



-[/oj^exp 

m,n 



— — {{am + x) 2 + (an + y) 2 ^j 



(15) 



Taking into account Eq. (15) and the following relation 

<U(x,y) >= -U Tr<7 2 /a 2 , 

we equate Eqs. (2) and (13) and obtain condition of 
the transition between the two vortex phases: 

= — - ■ — — ■ — — o ; AE 1 = Elat 



-(E™=-ooCX P (-^)) 



Ef, 



(16) 



This equation defines the boundary between the pinned 
lattice and the deformed triangular lattice. The phase 
diagram is presented in Fig. 3 in the Uo — a plane. Below 
the curve 1 the triangular lattice has the lower energy and 
upper - the pinned regular phase. The transition between 
these phases is discontinuous. It is clearly seen that the 
pinning efficiency depends not only on the amplitude of 
the potential but also on the size of the potential well. In 
fact it follows from Eq. (16) that the pinning efficiency 
is controlled by parameter 



Ua 



m— — oo 



exp 



The larger is the pinning site well the higher ampli- 
tude is required to stabilize the pinned phase. As we will 
see below, this phase diagram corresponds to the vor- 
tex lattice at all sub-matching fields not exceeding Hi 



/4j #3/8: ^1/3: I1 2/3 ), 



Hi/ 3 , H 2 



except H\ and Hi 



/2- 



when phase diagrams are different. 

At some specific matching fields the structure of vortex 
lattice can be more complicated. Below the first match- 
ing field, this happens at Hi and -ffi/2, when pinned lat- 
tice has a square symmetry (Fig. 1(d), (e)). In this case 
an intermediate phase appears with decreasing pinning 
strength, in which vortices in odd rows are depinned and 
are located in the middle between pinning sites, see Fig. 
4. Hence, only a half of the vortices is pinned. The exis- 
tence of this phase was predicted in Ref. |ll| at H = Hi 
for infinitely small potential well size. The symmetry of 
such lattice is close to the triangular one, and the en- 
hancement of the pinning energy is compensated by the 
reduction of the vortex interaction energy. We analyze 
now the phase diagram at H = Hi/ 2 , Hi taking into 
consideration the existence of this state. Similar to the 
previous cases the energy E int of the intermediate phase 
can be presented as 



Eint = E, 



1 + 1(U(0,0) + U(0,~)). (17) 



where E' t 

that E' tr is close to E t 
formed triangular lattice and the intermediate state is 
defined by the following equation: 



tr is the energy of the vortex interaction. Note 
tr ±o uuoc w uj tr . The boundary between the de- 



U = 2 (E' tr -E tr )x 




exp 

m— — oo 



(a/ay 



(18) 



The boundary between the intermediate phase and the 
square pinned lattice is given by: 




^ I a 2 (m + 0.5) 2S 
2^ exp I — | \ 



Z exp 



(19) 



where E sq is the energy of an ideal square vortex lattice 

{Elat = E sq ). 

The vortex phase diagram for Hi and Hi/ 2 and is pre- 
sented schematically in Fig. 5(a). This diagram is more 
complicated than that shown in Fig. 3. Below curve 1 
the deformed triangular lattice is more favorable energet- 
ically, in the region between curves 1 and 2 - the interme- 
diate state, and above curve 2 - the pinned square lattice. 
The transitions between all these phases are discontinu- 
ous. The values of E tr , E' tr , and E sq were calculated 
within the London theory at H = Hi , a = A for (A is the 
London penetration depth). The energy is measured in 
units of H 2 /8ir 7 where H c is the thermodynamic critical 
field. 

To analyze how the phase diagram depends on the type 
of the pinning potential, we also calculated the diagram 
for the parabolic pinning well: 



V(r) = -U a 



l-{ r - 
a 



if r < a; V(r) = 0, if r > a, 
(20) 



where r is the distance from the center of the well. In 
this case the following expressions define the boundaries 
between the triangular lattice and intermediate phase, 
intermediate phase and square lattice, triangular lattice 
and square lattice, respectively: 



U 



2 {E't r - E t; 
1 — 7TTT 



U a — 2 {E sq - E tr ) , 



(21) 



(22) 
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Uo = ^-^k. (23) 

The phase diagram for this potential is presented in Fig. 
5(b) at a = 0.2A for H = H\ There is also the interme- 
diate phase separating the deformed triangular and the 
pinned square lattices. The boundaries are represented 
by curves 1, 2, and 3. The main difference from the Gaus- 
sian interaction between the vortex and the pinning site 
is the existence of the triple point, when the size of the 
potential well is large. 

We can conclude that for any type of the pinning po- 
tential the phase diagram of the vortex lattice at H\ 
and H 1 / 2 in the plane pinning strength-pinning well size 
consists of three main regions: the deformed triangu- 
lar lattice (most vortices are situated not at the pinning 
sites), the ordered intermediate phase (half of vortices 
are pinned by PPA and half are interstitials) , and the 
pinned square lattice (all vortices are pinned by square 
PPA). However, details of the phase diagram depend on 
the exact form of the potential profile. For some types 
of the pinning potential the phase diagram has a triple 
point between these three phases, when the sizes of the 
potential wells are large. 

III. CRITICAL CURRENT 

In this subsection we study the critical current of the 
sample in the intermediate phase. We consider two situ- 
ations - when the external current is directed along the 
y- or the x-axis on Fig. 4(b). It is obvious that criti- 
cal currents are different in these cases in contrast to the 
pinned square and the deformed triangular lattices, when 
they are the same in the x- and y-directions. 

In fact there are two systems of the potential wells in 
the intermediate phase. The first one is formed by the 
pinning sites themselves and positions of these wells are 
fixed. The second one is formed due to the interaction 
with vortices trapped by pinning sites ("caging"). In 
absence of any current the second potential wells have 
minima at the interstitial positions, see Fig. 4(b). When 
the transport dc current is applied, the vortices trapped 
by pinning sites are displaced, and as a result, the po- 
sitions of minima of the potential wells, in which inter- 
stitial vortices are situated, will also change. Therefore, 
displacements of both interstitial and pinned vortices de- 
pend on each other and on the parameters of the pinning 
potential. 

Consider first the situation when the external current 
j flows in the y-direction and H = Hi. The Lorentz 
force acting on vortices is directed along the x-axis and 
is given by: 

h = (24) 

K 

Here and below the dimensionless variables are used: the 
distance r, the magnetic flux density h, and the energy 



density are measured in units of A, H c \/2, H^/8tt, re- 
spectively. In this notation $0 = 2it/k, where <I>o is the 
magnetic flux quanta. To calculate the energy of the sys- 
tem we use the usual London expression for the vortex 
field (h}: 

b(r) = -K (r), (25) 

K 

where Ko(r) is the modified Bessel function, r is the dis- 
tance from the vortex center. The magnetic energy of 
the vortex is given by [ ftpf : 

E m = —h{Q), (26) 

K 

where h(0) is the total magnetic field in the center of 
vortex. This field is equal to the sum of fields of each 
vortex. To describe the pinning potential we use the 
Gaussian potential defined by Eq. (15). We denote the 
displacements of interstitial and pinned vortices as xi 
and X2 ■ The energies Fi and F2 and of these vortices are 
given by: 




+IlX2- (28) 

We omitted all terms independent on X\ and x-i in these 
expressions. The vortex positions on x\ and Xi can be 
found from the force balance condition: 



These equations can be solved numerically. As a result 
we can find the critical current, i.e., the maximum value 
of the current, at which Eqs. (29) still have solutions. 
The critical current dependence on the amplitude of the 
vortex-pinning site interaction Uq is plotted in Fig. 6 at 
a = 0.1a, a = A. Below solid curve the vortex lattice is 
stable, above it becomes unstable and vortices start to 
move. 

The situation is more complicated, if the external cur- 
rent flows in the x-direction, see Fig. 4(b). For this case 
the system of equations similar to Eqs. (27)-(29) has 
been solved. Our results are presented in Fig. 7. We 
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found that the vortex depinning transition occurs in two 
steps. First, the positions of interstitial vortices become 
unstable at some finite value of current, and they jump 
to the neighboring pinning sites in the y-direction. The 
dependence of this current on Uq is shown in Fig. 7 at 
a = 0.1a, a = 1 (curve 1). The pinned square lattice 
remains stable up to the value of critical current (curve 
2) , which is more than two times higher than the current 
leading to the first instability. Thus, the phase diagram 
(Fig. 7) consists of three regions: phase with interstitial 
vortices, pinned square lattice, and moving vortices. The 
dependence of the critical current on Uq for this case is 
also presented in Fig. 6 (dot curve) for the comparison 
with the critical current in the y-direction (solid curve). 
The critical current in the first case is much higher than 
that in the second case. This is because for the destruc- 
tion of the stable vortex configuration in the first case we 
need to suppress the potential wells created by pinning 
sites. In the second case - only effective potential wells 
created by the repulsion with the pinned vortices. Such 
a high anisotropy of the critical current is a specific fea- 
ture of the intermediate phase, and it can be used for the 
experimental detection of this state. 

Finally, we note that according to our estimates, based 
on the results of Refs. the pinned phases are en- 

ergetically more favorable in experiments on the 

films with antidot arrays, since the pinning efficiency is 
very high in such systems. The triangular deformed lat- 
tice and intermediate phase can be studied in the films, 
where the pinning strength of the sites can be tuned, for 
instanse, by using the periodic arrays of blind holes. 

IV. CONCLUSIONS 

We have studied the stable vortex configurations in a 
regular array of weak pinning sites at matching fields not 
exceeding the first matching field H\. When the pinning 
efficiency is low enough it is energetically favorable to 
have a deformed triangular lattice, when it is sufficiently 
high a lattice of pinned vortices is formed. The phase di- 
agrams of the vortex lattices are calculated in the plane 
of two coordinates: the amplitude of vortex-pinning site 
interaction and the pinning potential well size. The vor- 
tex displacements in the deformed triangular lattice with 
respect to the positions in an ideal one are found within 
the elasticity theory. When the external field is equal to 
the first matching field or is twice lower, the existence of 
the intermediate regular phase is also energetically favor- 
able, in which half of the vortices is pinned by periodic 
pinning array and the vortex lattice is close to the tri- 
angular one. We analyzed the phase diagrams for this 
case taking into account different types of the pinning 
potential profile. We found that for some types of the 
potential the triple point can exist on the phase diagram, 
where deformed triangular lattice, pinned square lattice, 
and intermediate phase coexist. The critical current in 



the intermediate phase is calculated. This current turns 
out to be highly anisotropic. When current flows parallel 
to the pinned vortices rows the depinning occurs when 
the Lorentz force suppresses the effective potential wells, 
in which caged interstitial vortices are trapped. When 
current flows perpendicular to the pinned vortices rows, 
first these effective "caging" potential wells are destroyed, 
and interstitial vortices jump to the neighboring vacant 
pinning sites. The phase diagrams for the stable vortex 
lattice in the plane of applied current-pinning potential 
amplitude are also obtained. 
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FIGURE CAPTIONS 



FIG. 1. Vortex patterns corresponding to different matching fields: (a) - 1/8H\ , (b) - 1/AHi, (c) - 1/3-ffi, (d) - 
l/2Hi, (c) - Hi. Open circles show positions of the pinning sites, black dots - vortices. 

FIG. 2. The structure of the triangular vortex lattice, deformed by the square pinning array of Gaussian potential 
wells, defined by Eq. (20), at H = H 1 , c n S /U = 0.2, a = 0.1a. 

FIG. 3. Phase diagram in the plane of the vortex-pinning site interaction amplitude and the pinning potential 
length-scale at sub-matching fields not exceeding Hi, except and (-Hi/s, Hi/ 4, H 3 / & , -H1/3, ^2/3)- Below curve 1 the 
deformed triangular lattice has the lowest energy, above -pinned phase. 

FIG. 4. The structure of the intermediate phase at #1/2 (a) and Hi (b). Black dots and open circles denote the 
positions of the vortices and the pinning sites. Dashed lines show the symmetry of the vortex lattice. 

FIG. 5. Phase diagram in the plane of vortex-pinning site interaction amplitude Uo and the pinning potential 
length-scale a at H = Hi. Fig. 5(a) corresponds to the Gaussian potential of the vortex-pinning site interaction, (b) 
- to the parabolic one. The similar phase diagrams are also valid for H = Hi/ 2 . 

FIG. 6. Phase diagram of the intermediate phase in the plane: applied current - amplitude of the vortex-pinning 
site interaction at H = Hi, a = 1, a = 0.1a. The current flows along the y-direction in Fig. 4(b). Solid curve (critical 
current) shows the boundary between the stable lattice with the interstitial vortices and the moving vortices. Also 
shown for the comparison (dot curve) is the critical current in the ^-directions in Fig. 4(b). 

FIG. 7. Phase diagram of the intermediate phase in the plane: applied current - amplitude of the vortex-pinning 
site interaction H = Hi, a = 1, o = 0.1a. The current flows along the ^-direction. Curve 1 corresponds to the 
transition from the lattice with interstitial vortices to the pinned square lattice, curve 2 - to the depinning transition. 
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Fig. 1(a) 
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Fig. 1(b) 



8 



o 


® 


o 


o 


® 


o 












(c) 


o 


o 


® 


o 


o 


® 


® 


o 


o 


® 


o 


O ' 


o 


® 


o 


o 


® 


o . 


o 


o 


® 


o 


o 


® 


® 


o 


o 


® 


o 


o 



9 



Fig. 1(c) 
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Fig. 1(d) Fig. 1(c) 
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